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Abstract 

The neutral Higgs sector of the next-to-minimal supersymmetric standard model (NMSSM) 
with explicit CP violation is investigated at the 1-loop level, using the effective potential 
method; not only the loops involving the third generation of quarks and scalar quarks, but 
also the loops involving W boson, charged Higgs boson, and chargino are taken into account. 
It is found that for some parameter values of the NMSSM the contributions from the W 
boson, charged Higgs boson, and chargino loops may modify the masses of the neutral Higgs 
bosons and the mixings among them significantly, depending on the CP phase. In e^e~ 
collisions, the prospects for discovering neutral Higgs bosons are investigated within the con- 
text of the NMSSM with explicit CP violation when the dominant component of the lightest 
neutral Higgs boson is the Higgs singlet field of the NMSSM. 
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I. INTRODUCTION 



In the standard model, the phenomenology of CP violation is described in terms of a complex 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix for quarks [1]. The source of 

the complex phase is not clarified in the standard model. If there are two Higgs doublets in the 
model, the complex phase might originate from the Higgs sector. As the supersymmctrically 
extended versions of the standard model require at least two Higgs doublets [2] , to give separate 
masses for the up- type quarks and the down-type quarks, they are regarded as natural candidates 
for explaining the origin of CP violation. This is all the more because supersymmetry is broken 
in nature at the electroweak scale. One would obtain no CP phase in the Higgs sector if 
supersymmetry is intact. 

In supersymmetric standard models, the CP violation can occur in the neutral Higgs sector 
through the mixing between the scalar and pseudoscalar Higgs bosons. In principle, it can 
be accomplished either spontaneously or explicitly. Either way, the phenomenology of the CP 
violation in the supersymmetric standard models is considered as particularly useful probes of 
new physics at the high energy scale because CP violation only appears through the CKM matrix 
in the standard model. Recently, the attempts to explain the source of the CP violation within 
the context of supersymmetry have been investigated by many authors [3-5]. 

In the minimal supersymmetric standard model (MSSM), many authors have investigated 
if it is possible to violate the CP symmetry either spontaneously [3] or explicitly [4,5]. It is 
found that the tree level Higgs potential of the MSSM is unable to break the CP symmetry in 
either way. Further investigations at the 1-loop level of the MSSM show that the scenario of 
spontaneous CP violation is experimentally ruled out by the Higgs search at the CERN e^e~ 
collider LEP, because it leads to a very light neutral Higgs boson [3]. However, at the 1-loop level, 
the scenario of explicit CP violation in the MSSM is found to be viable through the radiatively 
corrected Higgs potential. Subsequently, a number of investigations have been devoted to the 
effects of explicit CP violation in the MSSM at the 1-loop level to the charged Higgs boson 
as well as to the neutral ones by considering radiative corrections due to the third generation 
quark and scalar quark contributions [4]. Recently the MSSM searches for explicit CP violation 
were extended to radiative corrections due to the W boson, the charged Higgs boson, and to 
the chargino contribution [5]. 

Meanwhile, the CP symmetry may also be broken if the Higgs sector is extended by intro- 
ducing at least one Higgs singlet. The supersymmetric standard models with additional neutral 
Higgs singlets can evade some of the strong constraints of the MSSM. These models may be 
classified according to whether they have a discrete symmetry or not. It is known that for those 
models without a discrete symmetry, spontaneous CP violation is possible at the tree level 
[6] and up to the 1-loop level including radiative corrections due to the top and scalar top quark 
loops [7]. On the other hand, those models with discrete symmetry can cause a domain problem 
associated with the existence of degenerate vacua after the electroweak phase transition [8] . The 
possible domain wall problem is assumed to be solved by adding nonrenormalizable interactions 
which break the symmetry without spoiling the quantum stability [9]. Among those with a 
discrete Z^ symmetry on its Higgs potential, the simplest model, with just one Higgs singlet, 
is usually called as the next-to- minimal supersymmetric standard model (NMSSM) [10]. The 
NMSSM is basically different from the MSSM in that it has one extra Higgs singlet. 

For spontaneous CP violation, it is well known that the NMSSM, which has Z3 symmetry 
in its Higgs potential cannot produce spontaneous CP violation at the tree level, because of 
vacuum stability [11]. This behavior does not change by taking into account the contributions 
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of the degenerate scalar top quark loops [12]. Spontaneous CP violation can only occur if the 
scalar top quark masses at the 1-loop effective potential are nondegenerate [13]. However, at the 
tree level, unlike the MSSM, the scenario of explicit CP violation is possible for the NMSSM. 
By assuming the degeneracy of the top scalar quark masses in the Higgs sector of the NMSSM, 
it is found that large explicit CP violation may be realized as the vacuum expectation value 
(VEV) of the neutral Higgs singlet approaches the electroweak scale [14]. 

In a previous article [15], we have analyzed the phenomenology of explicit CP violation in the 
Higgs sector of the NMSSM at the 1-loop level. Using the effective potential method, we have 
included the radiative corrections due to quarks and scalar quarks of the third generation. We 
have found that there arc parameter regions of the NMSSM at the 1-loop level where the lightest 
neutral Higgs boson may even be massless, without being detected at LEP2. The VEV of the 
neutral Higgs singlet has been shown to be no smaller than 16 GeV for any parameter values of 
the NMSSM with explicit CP violation. This value of the lower limit is found to increase up to 
about 45 GeV as the ratio (tan /?) of the VEVs of the two Higgs doublets decreases to smaller 
values (~ 2). The discovery limit of the Higgs boson search at LEP2 is found to cover about 
one half of the kinematically allowed part of the whole parameter space of the NMSSM, and the 
portion is roughly stable against the CP phase. 

In this article, we continue our analysis more exhaustively by including the contributions from 
the loops involving the W boson, charged Higgs boson, and chargino, in addition to those from 
the loops involving quarks and scalar quarks of the third generation. It is found that the inclusion 
of the W boson, charged Higgs boson, and chargino loops yields significant modifications to the 
results on the neutral Higgs sector obtained without them, depending on the parameter values 
of the model. At the 1-loop level two CP phases arise from each scalar quark mass matrix 
for the third generation. We assume that at the tree level the CP phase is equal to two CP 
ones arising from the contribution of radiative corrections due to the third generation quark 
and scalar quarks. One CP phase appears from the 1-loop effective potential coming from the 
chargino sector. We assume that the CP phase arising from the chargino sector does not equal 
the other ones at the 1-loop level. We investigate the size of the contribution of the W boson, the 
charged Higgs boson, and the charginos to the neutral Higgs boson mass at the 1-loop level. We 
then investigate the phenomenological consequences for the neutral Higgs sector of the NMSSM 
with explicit CP violation at future e+e" collisions. 

II. EFFECTIVE POTENTIAL 

The additional Higgs singlet superfield in the NMSSM can avoid the so-called yU-parametcr 
problem in the MSSM, where the dimensional ii parameter of the MSSM is introduced by hand. 
In the NMSSM, the corresponding quantity is generated dynamically. The additional Higgs 
singlet superfiled, N, has a vanishing hypercharge and neutral. The neutrality of N enlarges 
the spectrum of neutral particles of the NMSSM with an additional scalar Higgs boson, an 
additional pseudoscalar Higgs boson, and an additional neutralino, while the structure of the 
charged sector of the NMSSM is essentially the same as that of the MSSM. 

Ignoring all quark and lepton Yukawa couplings except for those of the third generation 
quarks, the superpotential involving the superfields Hf = [H^,H^), HJ = (iJ^,iJ2), and N in 
the NMSSM with discrete symmetry is given by 

W = hbHleQb% - htH^eQt% + \H'(eH2N - ^N^ , 
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where ht and /i5 arc Yukawa couphng constants of the quarks for the third generation, and the 
parameters A and k arc dimcnsionless couphng constants. The 2x2 antisymmetric matrix e is 
defined as ei2 = —621 = 1. The SU(2) doublet superfield = {tL,bL) contains the left-handed 
quark for the third generation, and the SU(2) singlet superfields f^^ and 6^ are the charge 
conjugates of the right-handed quarks for the third generation. At the tree level, a massless 
pseudoscalar Higgs boson would appear if the global U(l) Peccei-Quinn symmetry would not be 
broken explicitly by the cubic term in the Higgs singlet in the superpotential. 

In the NMSSM one breaks the supersymmctry by introducing explicitly soft breaking terms 
to the Lagrangian. The relevant soft breaking terms for the above superpotential can be written 
as 

Aoft = -ml,J/fi|2-ml,Ji?2|'-m?,|iV|2-m||Qp-mi.|t^p-m||6^|2 
+ {-hAbHfeQb% + htAtH^eQi% + H.c.) 

+ (^XAxHfeH2N + ^AkN^ + H.c.) + QmsA^A^ + ^MiAiAi + H.c.) , (1) 

where m^_^, rnj^^, rnj^, niQ, m^, and are the soft breaking masses; Ax, A^, At, and Af, are 
the trihnear soft SUSY breaking parameters with mass dimension; Ag (a = 1,2,3) and Ai are 
the gauginos of the SU(2) and U(l) gauge groups; and M2 and AIi arc their mass parameters, 
respectively. The full tree level Higgs potential can then be obtained by collecting the relevant 
F terms and D terms in the superpotential, to add to the Higgs part of the above soft breaking 
terms. 

The condition for the electroweak symmetry breaking is obtained by minimizing the Higgs 
potential V = V^' with respect to the VEVs of the three Higgs fields. We assume that the 
VEVs of three neutral Higgs fields H^, H2, and N are, respectively, vi, V2, and x, all being real 
and positive, so that there is no spontaneous CP violation in the Higgs sector in the NMSSM 
[11-13]. The ratio of the two VEVs is denoted by tan /3 = V2/V1, giving the Z boson mass as m| 
= {gf + 52)^^/2, and the W boson mass as = g^v"^ /2, with v = y'v'l+v'l = 175 GeV. The 
soft SUSY breaking masses m^^, vn^^, and can be eliminated by imposing the minimum 
conditions on the Higgs potential. 

The explicit CP violation in the Higgs sector of the NMSSM might be activated if there 
are complex phases in the Higgs potential. From the tree level Higgs potential of the NMSSM, 
complex phases might occur by assuming that A, Ax, k, and A^ are complex numbers. Among 
them, XAx and kA^ can be adjusted to be real and positive by redefining the overall phases of 
the superfields Hi, H2, and N. Then, we are left essentially with one phase: It may be chosen 
to be the one in Xk* = A/ce*'^" , which allows for the explicit CP violation at the tree level. This 
phase appears in the expression for the mass of the charged Higgs boson [15]: 

2x 

= rriw - X^v^ + X{Ax + kx cos ^0) . ■ (2) 

Sin 

Note that although the phase ^0 is present in the potential, it does not show up in the vacuum. 
Thus, one need not impose a minimum condition on the Higgs potential with respect to it. 
Consequently, the tree level Higgs potential contains ^q, tan/3. A, Ax, k, A}^, and x as free 
parameters. 

Now, at the 1-loop level, in the NMSSM with explicit CP violation, the radiative corrections 
to the tree level Higgs potential come primarily from the quark and scalar quark loops of the 
third generation: Generally, the leading radiative corrections come from the top and scalar top 
quark loops. For large tan/3, the contributions of the bottom and scalar bottom quark loops 
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can be significantly large. Nevertheless, in order to carry out a comprehensive study on the 
phenomenology of the explicit CP violation at the 1-loop level in the NMSSM, one should also 
include the contributions of the W boson, the charged Higgs boson, and the charginos loops. 
Let us decompose radiative corrections to the 1-loop effective potential into three parts: a part 
for the radiative corrections from the top and scalar top quark loops, another one for those from 
the bottom and scalar bottom quark loops, and the third one for the radiative corrections from 
the W boson, the charged Higgs boson, and the charginos. Thus, the radiative corrections to 
the 1-loop effective Higgs potential is written as [16] 



= V* + V'' + , (3) 



with 



V' = 




(4) 

where Ai denote the mass matrices for the particles and superparticles with the field dependence, 
and A is the renormalization scale in the modified minimal substraction (MS) scheme. The full 
1-loop level Higgs potential will be given by adding to the tree level Higgs potential. 

Among the mass matrices, there is a 4 x 4 Hermitian mass matrix for the scalar quarks 
of the third generation, which may be broken up into four 2x2 block submatriccs. After the 
electroweak symmetry breaking, only two diagonal block submatrices become dependent on the 
neutral Higgs fields and thus have nonzero eigenvalues. By diagonalizing each diagonal block 
submatrix, one can obtain the 1-loop masses of the scalar top quarks and the scalar bottom 
quarks separately. From the complexity of the relevant parameters, a phase appears in the 
expressions for the masses of the scalar top quarks and another one for the masses of the scalar 
bottom quarks. These two phases may also take part in the explicit CP violation at the 1-loop 
level. In general, the two phases are different. However, for simplicity, we assume that Aj, = Af, 
while keeping them complex. Then, the two phases become equal. 

The masses of the scalar top quarks m?^ and m? arc obtained as 



1-1 



rrii 



+ hi'^Q + "^t) T \l\{'mQ - rn^Y + m1{Af + \^x^ cot^ /3 + 2AtXx cot /3 cos cf)) , (5) 



and those of the scalar bottom quarks m? and m? as 

01 02 



ml + \{mQ + rriT) T \I\{^q - + "^l^M + ^^^^ ^an^ (3 + 2At\xl&nf3 cos (j)) , (6) 

where the contribution from the D terms are not included. Thus the CP phase in the bottom 
scalar quark masses is equal to that of the top scalar ones: (f) is the relative phase between At = A^ 
and A, allowing for the explicit CP violation in the Higgs sector through radiative corrections. 
The phase 4> would be absent as long as both the scalar top quarks and scalar bottom quarks 
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are degenerate in masses. In the scalar quark masses (f) has no CP- violating effects on the Higgs 
sector, since it could be replaced by some modification At anyhow. However, unless vanishes, 
explicit CP violation occurs in the Higgs couplings through radiative corrections. We assume 
that (/) is equal to (?!)o, in the mass of the charged Higgs boson, for simplicity. Note that the top 
quark mass is given by rrif = h^v"^ sin^ /3, and the bottom quark mass as = h^v'^ cos^ /3. 
The square masses of charginos are given by the eigenvalues of the square-mass matrix 

2 _ / |M2p+52|//,|2 g,{M2H^,*-X*HlN* 
as 

^li, 2 - H^l + ^^^^ + 2"^^} T [i{M| - A\2}2 + ml,{Mi - X^x^} cos 2(3 

+ 2m^{(M2sin/3)2 + {Xxcospf - M2Axsin2/3cos0c}]^''^ , (7) 

where we neglect the quartic term of g2 since its order is equal to that of the D terms in the 
scalar quark sector. Notice that there is a complex phase in the chargino masses, (j)c, arising 
from the relative phase between the SU(2) gaugino mass M2 and A. We assume that 0c is not 
equal to 0. The total Higgs potential at the 1-loop level is written as V = V'^ + V^. Here, 
contains (pc, tuq, itit, At, and M2 as free parameters. 

The Higgs sector of the NMSSM has ten real degrees of freedom coming from two Higgs 
doublets and one Higgs singlet. After the electroweak symmetry breakdown takes place, three 
of them correspond to a neutral Goldstone boson and a pair of charged Goldstone bosons, and 
the other seven correspond to five neutral Higgs bosons and a pair of charged Higgs bosons. 
Transforming to a unitary gauge, one may express two Higgs doublets and one Higgs singlet as 



H 



H, 
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vi + Si + i sin (3 A 
sin/3C+* 

cos/?C+ \ (8) 
V2 + S2 + icos/3A ) ' 
AT = { x + X + iY ) , 

where ^i, ^2, A, X, and Y are neutral fields and C+ is the charged field. If CP symmetry is 
conserved in the Higgs potential. Si, S2, and X would become the CP even scalar Higgs fields, 
whereas A and Y would become the CP odd pscudoscalar Higgs fields. The three orthogonal 
combinations among A and yield one neutral Goldstone boson and a pair of charged Gold- 
stone bosons; they will eventually be absorbed into the longitudinal component of Z and W 
gauge bosons, respectively. 



III. NEUTRAL HIGGS BOSON MASSES 

In the NMSSM with explicit CP violation the neutral Higgs boson mass matrix is given as a 

5x5 matrix that is obtained by the second derivatives of the Higgs potential with respect to 
the five neutral Higgs fields. On the basis of {Si,S2,A,X ,Y) the elements of the symmetric mass 
matrix for the neutral Higgs boson can be expressed as 

Mij = + 6Mlj + dM^j + 5Mf. (i, j = 1 to 5) , (9) 

where we note that Mf^ comes from the tree level Higgs potential, SMj^ from the contribution 
of the top and scalar top quark loops, 5M^, from the contribution of the bottom and scalar 



6 



bottom quark loops, and SM^j from the loops of the W boson, the charged Higgs boson, and 
the charginos. Thus, MP. is obtained from 5Mf- from F*, 6MP- from and SM^^ from V^. 

The complicated expressions for Mf-^ ^^Ij^ 6M^- are given elsewhere [15]. We concen- 
trate on SMf-. We calculate each elements of it. The results read 

5M^^ = ^,22 T 2 2 A2 + W^^2Ax tan /3 cos (/)c/ m? , 

m^cos=^/3A^,log(m|^/m|J m^cos^/3 /m|^rn|+ \ 

2 _ ^2 \ + o^2„2 \ A A 



27r2-y2 C^2 _ j^2 \ StT^v"^ . ,.„, 

^ X2 XI ^ \ XI X2 

t2 Q ( / rm'^ 



r^2,2 (m|^-m|J2 +^M2Aa;cot/jcos<^J(m|^, m|j 
sin2 /? A;^2 log(m|^/m|J ^ sin^ /3 ^^^^2,+ ^ 



27r2ti2 (tt^I _ j^2 ^ Stt^w^ I jt^I ^4 

A,2 A,l ^ * ^ " 

+ ^f^(2A'^' - m^,) [ml, log f - m 



COS^ /3 / , 2 2 2 N f 2 1 / "^r 

^^^(2A ^; - m,^) log j - 

M| A^x^ sin^ 0c ,2 2 \ mp^M2AxcoS( 



^ X2 XI' 



X2 'Xi' 

_ / m|^m| \ A3.T A^log(m|^/m|J 

167r2 ^[ I 87r2 (m|^-m|^^ 



X2> 



>^ Axsm2P [ 2 \ _^2 

xj^^x m^M^X^sm^2(5s\v? ct)c ,2 2 ^ 

&5 47r2 (m|2-m|J2 'i:!' '<:2^ 

+ ^--^M2Asin2/3cos(;^c/("^|l, w-lJ 

Asin2/3^^ ^ ^, ^^ f 2 
327r2 



(Aa + Sfcxcos 0) |m^+ log (-^^ - "^c+ 



^^?2 = -5^^^^^^K.-|2)-5^M2Axcos0jK,,m|,] 

_ sin2/3(A;^, + A^ 

87r2u2 (ml^ 

sin 2/3 / m^m^+ \ A2m^sin2/3 f ml 
+ Zoo log — T— ^ SZ9 log ' 
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mfy M2XXC0S PsincpcAj^^ 



27r2u2 



i^x2 - "^i ) 



XI ^ 



M2Aa:cos/3sin^c/("i~., m~ ) 



+ 



m 



w 



M2X.r cos ;i sin (;6c 




A cos /SA-j^j 



Svr^t; (ml 



X2 



+ 



m 



w 



m^^X^x cos P 



XcosP{\x — M2 tan P cos (f>c)f{m'~.^, m-^) 



4. 



log 



A cos /3(AxAjj.j + Ajj.) 



A4 



log 



XI ^ 



M2A cos /3 sin 2/3 sin <;i!)c A 



47r^'y 



XI' 



^2) 



+ 



M2Asin/3sin(?!)c/(m|j, m^J 



rriy^ M2 A cos /3 sin 2/3 sin 



XI ^ 




M2Axsin/3sin(/>cA^2 2 



X2' 



m 



w 



+ 



M2XX sin P sin (pcfim-^^, m\^) 
log 



47^2^,2 

M2Aa; sin /3 sin 



X2 




Asin^Ap^jAj^ 2 



XI' 



X2^ 



XI ^ 



+ 



m 



w 



2 

w 



X sin /3{Xx — M2 cot /? cos 4>c) f {m-.^ , m-^) 



4. 



m^rX^x sin (3 



log 



Asin/3(AxA;^2 + Ax) 



A4 



m|j^ M2Asin/3sin2/3sin(/)cA^2 ^_2 



K2-"^li) 



log 



47r2'y 



K2-^li)' 



+ 



w 



M2Acos/3sin^c/("^|i, J^i-J 



47r2u 

M2 A sin /3 sin 2/3 sin 
rriy^ M2A2xsin0cA 



■log 



m 



X2 



m 



XI 



XI > 



m 



w 



47r2? 



M2Asin(/)c/(m? , m|j + 



mL M2A^a;2 sin 0c 



M| A^x sin 2/3 sin^ i^c / 2 2 



X2 



XI ' X2 4^2^ ^j^2^ _ j^2 ^ 




^-"^^45 - g^2 (^^2 _ „.2 N2 9\^xi, ^X2) 

myi^ M2\^x sin 2/3 sin < 



87r2 ("^la-^^li) 




where we define 



A;^j = Af| + A^x^ - 2M2 Ax tan (3 cos 0c , 

A^2 = Af| + A^x^ - 2M2 Ax cot (3 cos 4>c , (11) 
= A^x^ - M|Ax + 2m^Ax - 2m^M2 sin 2/3 cos 4>c , 



and 



2\ 



/(mi, 1712 



1 



— mf) 



2 2 
2 1 2 1 ^2 

mi log - "^2 log ^ 



+ 1 



mo + m? , 



5(m?,mi) = ^^^l^log^ + 2. (12) 

At the tree level, in M^?-, the scalar-pseudoscalar mixing between two Higgs doublets would not 
occur from the Higgs potential since both M^g and M23 are zero. Meanwhile none of the other 
elements for the scalar-pseudoscalar mixing in are zero. Assuming the nondegeneracy of the 
scalar quark masses, the radiative corrections due to the quark and scalar quark for the third 
generation yield the scalar-pseudoscalar mixing between two Higgs doublets. Note that at the 
1-loop level, on the contrary, every element of 6M^j is nonzero; the magnitudes of these elements 
are proportion to sin (pe- 
lt is impossible to obtain analytic expressions for the eigenvalues of Adij as functions of the 
various parameters. We carry out the job of obtaining the eigenvalues numerically. The masses 
of the neutral Higgs bosons are then sorted such that m^. ^ nih for i < j, such that m/^ is the 
mass of the lightest neutral Higgs boson. For the numerical calculations, wc set the ranges of the 
relevant parameters in the Higgs potential. The renormalization scale in the effective potential 
is taken to be 500 GeV. We fix rrit = 175 GeV for the top quark mass [17], mi, = 4 GeV for 
the bottom quark mass, mw = 80.4 GeV for the charged weak gauge boson mass, and mz = 
91.1 GeV for the neutral weak gauge boson mass. The upper bounds on A and k are estimated 
as 0.87 and 0.63, respectively, by the renormalization group analysis of the NMSSM [18]. Thus, 
we set their ranges as < A ^ 0.87 and < k ^ 0.63. The ranges for other parameters are set 
as follows: 2 ^ tan/3 ^ 40, < ^a, ^fc, x, mq, mx (= m^) ^ 1000 GeV, < At = Ab ^ 2000 
GeV, and < M2 ^ 500 GeV. Note that there are additional parameters at the 1-loop level, 
arising from V^: the relative phase {(j)c) between complex M2 and A, as well as <p. We assume 
that they may vary from to tt. We further assume that the lighter scalar quark masses of the 
third generation are larger than the top quark mass. 

Then, we impose phenomenological constraints on the values of the parameters. Prom the 
negative results at the LEP2 experiments on the neutral Higgs boson, research has already 
suggested some constraints on the parameter space of the NMSSM with explicit CP violation 
[15]. At the center of mass energy of LEP2, -s/s = 200 GeV, the dominant production processes 
for the neutral Higgs boson are (i) the Higgsstrahlung process, e+e~ Z* ^ Zhi {i = 1 to 5) 
and (ii) the Higgs-pair production process, e'^e~ Z* ^ hihj = 1 to 5, i ^ j). The total 



9 



cross section for producing any one of the neutral Higgs bosons in e+e^ collisions is then defined 
by the sum of the two processes: at = o"i + (Jij/2 for i ^ j. Taking the discovery limit 
of LEP2 as 0.1 pb for ^/s = 200 GeV, we exclude the parameter values if they produce at ^ 0.1 
pb. 

We introduce a useful expression in order to evaluate the effects of explicit CP violation on 
the neutral Higgs boson masses as in the case of the MSSM with CP violation in the Higgs 
potential [19]. It is the dimensionless parameter given by 

p = B'Ol.OlOlOlOl , (13) 

where Oij = 1 to 5) are the elements of the orthogonal transformation matrix that diagonal- 
izes the neutral Higgs boson mass matrix. The range of p is from to 1 since the elements of the 
transformation matrix satisfy the orthogonality condition of Yl^=i Oji = 1- If /? = 0, there is no 
explicit CP violation in the Higgs sector of the NMSSM. On the other hand, if p = 1, CP symme- 
try is maximally violated through the mixing between the scalar and pseudoscalar Higgs bosons. 
The maximal CP violation that leads to p = 1 takes place when Of^ = O21 = O31 = O41 = O51 
= 1/5. The elements Oij determine the couplings of the physical neutral Higgs bosons to the 
other states in the NMSSM. Elsewhere, the chargino sector contributions, i.e., the contributions 
from the loops of the W boson, the charged Higgs boson, and the charginos, to the neutral Higgs 
boson masses as well as to the mixings among them have been analyzed within the context of 
the MSSM with expHcit CP violation [5]. 

Let mfej denote the full 1-loop mass of the lightest neutral Higgs boson, including the contri- 
butions from the loops of the W boson, the charged Higgs boson, and the charginos, as well as 
other loops, and denote the corresponding mass including all the loop contributions except 
the loops of the W boson, the charged Higgs boson, and the charginos. Therefore, m^^ is not 
the tree level mass but the 1-loop mass with the radiative corrections due to the top quark and 
scalar top quark loops as well as the bottom quark and scalar bottom quark loops. Only the 
chargino sector contributions are absent in it. Then, the mass difference between m/j^ and m^^, 
defined as 

would evidently express the amount of the contributions from Mf-, and hence V^, to m/j^. 

We scan the parameter space of the NMSSM using the Monte Carlo method and select points 
that are consistent with the above constraints. The selected points are physically allowed in the 
sense that at those points the Higgs boson masses are positive and at is smaller than 0.1 pb. 
Those points consist of the allowed region in the parameter space of the NMSSM. We search in 
the allowed region for several points where m^^ is large and/or p is large, i.e., where the effect 
of the chargino sector contributions to the smallest mass of the neutral Higgs bosons and/or 
the mixing among them is large. Figure 1 shows three such points with large m^^ and/or p 
as illustration. In Fig. 1, we plot m^^, nih^, \m^^\/mh^, and p, against the CP phase (j)c- For 
l/m/jj and p, we express them in %. The parameter values for Fig. 1(a) are cp = tt/2, tan/3 
= 3' A = A; = 0.5, = 800 GeV, = At = 1000 GeV, x (= my) = 500 GeV, mg = 900 
GeV, and M2 = 300 GeV, while the parameter values for Fig. 1(b) are cf) = 7r/2, tan/3 = 10, A 
= 0.03, Ax = 60 GeV, k = 0.4, A^ = 200 GeV, x = 500 GeV, tuq = mr = At = 600 GeV, M2 
= 100 GeV, and for Fig. 1(c) are (j) = tt/2, tan/3 = 10, A = 0.7, = 400 GeV, k = 0.5, Ak = 
20 GeV, X = 50 GeV, mq = 500 GeV, tut = At = 1000 GeV, M2 = 100 GeV. 

For the parameter values of Fig. 1(a), it is found that all of the neutral Higgs bosons m/j. 
(i = 1 to 5) as well as at are roughly stable against the variation of (/>£. They are numerically 
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obtained as m/j^ = 25, = 127, nihs = 786, nih^ = 819, nih^^ = 902, m(j+ = 815 GeV, and 
cr^ = 0.5 fb for = 200 GeV. One can see in Fig. 1(a) that the curve for m^^ = 42 GeV 
is flat. This is because m^^ contains neither (pc nor M2. Note that the effective potential for 
the contributions of the quark and scalar quark loop for the third generation is independent of 
both (pc and M2. In Fig. 1(a), rn^_ is always negative for the whole range of (pc Skud slightly 
decreasing. Thus, nih^ (GcV) is roughly stable but decreases quite slightly as 0c increases from 

to TT. Hence, the slight decrease in \m^^\/mhi for larger (pc- It is observed that p increases 
as 4>c goes from to tt. This behavior of p arises from destructive interferences between the 
chargino sector contributions and the other ones. 

Figure 1(b) displays a relatively large m^jj. We have m/^^ = 94, = 98, mh.j^ = 128, 

m/j4 = 346, m/jg = 346, m(j+ = 124 GeV, and at = 36 fb. For the parameter values of Fig. 
1(b), the neutral Higgs boson masses as well as at are almost stable against (pc, similar to Fig. 
1(a). One sees that and m/jj are nearly equal, and large as compared to the case of Fig. 
1(a). Thus, m^^ for the parameter values of Fig. 1(b) is very small. It is positive for smaller 
(pc and negative for larger (pc- There is a crossover at about (pc ~ 27r/3. On the other hand, 
p shows interesting behavior. It varies quite abruptly; it drops from 100% at 0c = down to 
zero at about <pc ~ O.SStt, and then increases as large as 80%. However, it is difficult to analyze 
why p varies so widely mainly because of the nature of the Higgs sector of the NMSSM. The 
complicated form of the 5x5 mass matrix of the neutral Higgs bosons prohibits us to trace 
analytically down the behavior of p. 

For the parameter values of Fig. 1(c), m^^ is always positive but decreasing from 13 GeV to 
5 GcV as (pc goes from to tt. Both m^^ and nih^ in Fig. 1(c) arc comparable in size to the case 
of Fig. 1(a). Approximately, we have m/i^ = 38, rrih^ = 52, m/jg = 124, m/j^ = 392, m/jg = 403, 
777,(7+ = 364 GeV, and at = 10 fb in Fig. 1(c). For m°^, we have about 25 GeV. The behavior 
of p for the parameter values of Fig. 1(c) is as varied as that of the case of Fig. 1(b). We note 
that the three sets of parameter values used in Fig. 1 arc consistent with the LEP2 data, as 
indicated by the size of the total cross sections for Higgs production: None of at is larger than 
0.1 pb, the discovery limit of LEP2. 

Now, let us study the behavior of m^^ in detail. It comes from the 1-loop effective potential 
V^, which includes the effects of the chargino mass splitting and possesses the CP phase (pc, 
which in turn depends on the SU(2) gaugino mass M2. It represents the contributions of the 
loops of the W boson, charged Higgs boson, and the charginos. Thus, the dependence of m^^ 
on (pc and M2 would give us a measure of the effect of explicit CP violation in the NMSSM. 

Figure 2 shows the dependence of on them, plotted against the ratio of the two chargino 
masses m,j^.^/m,j^^. In Fig. 2 there are three sets of curves, and each set consists of three curves 
corresponding to three different values of ^c^ 0c = (solid curve), 7r/2 (dashed curve), and tt 
(dotted curve). The extent of each curve corresponds to the variation of M2: We vary M2 from 

1 to 500 GeV. For the rest of relevant parameters, we set their values as: (p = 7r/2, tan/3 = 3, 
A = A; = 0.5, Ax = 800 GcV, Ak = At = 1000 GeV, x (= mr) = 500 GeV, and mq = 900 GeV 
for the set of curves of feature (Al); (p = 7r/2, tan/3 = 10, A = 0.03, Ax = 60 GeV, k = 0.4, 
Ak = 200 GeV, x = 500 GeV, and mg = uit = At = 600 GeV for the set of curves of feature 
(A2), and (p = 7r/2, tan^ = 10, A = 0.7, = 400 GeV, k = 0.5, Ak = 20 GeV, x = 50 GeV, 
niQ = 500 GeV, and rriT = At = 1000 GeV for the set of curves of feature (A3). In short, the 
parameter values for the three sets of curves of features (Al), (A2), and (A3) are the same as 
Figs. 1(a), 1(b), and 1(c), respectively, except for 0c and M2. 

When both 0c and M2 approach zero as the charginos becomes degenerate in mass, Tn^^/m^^ 
approaches 1. Still in this case, if 0c = 0, CP violation in the Higgs sector is possible from 
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yO _^ yt _^ yb gjj^^jg ^ fixed as tt/2. For the parameter values of feature (Al), m^^ is shown 
to be relatively independent of the variation in (^c; and found to increase as increases. 
Nevertheless, m^^ is always negative from —24 to about —16 GeV for the whole range of M2. 
The ratio of the chargino masses decreases from about 15 to 2 as M2 increases from 1 to 500 
GeV for the parameter values of feature (Al). On the contrary, for the parameter values of both 
features (A2) and (A3), the ratio of the chargino masses increases as M2 increases from 1 to 500 
GeV: For those of feature (A2), the ratio of the chargino masses increases from 1 to as large as 
30, and for those of feature (A3), from 1 to 15. 

The variation in m^^ for the parameter values of feature (A2) is rather small: ±2 GeV for 
the whole ranges of < M2 500 GcV and of ^ ^ vr. In contrast, the variation in 
for the parameter values of feature (A3) is very large: it can be as large as 16 GeV. Therefore, 
the overall view of Fig. 2 implies that m^^ is not negligible at all. Our analysis shows that for 
some particular values of the relevant parameters m^^ may have any value between —24 and 16 
GeV. Thus, it may play a quite crucial role in the analysis of the explicit CP violation scenario 
of the NMSSM. 



IV. NEUTRAL HIGGS PRODUCTION IN e+e COLLISIONS 

We know that the unsuccessful result at LEP2 for the Higgs search does not exclude the possibil- 
ity of the existence of a massless neutral Higgs boson in the NMSSM with explicit CP violation. 
In this section we investigate the detectability of a neutral Higgs boson in e'^e~ collisions with 
much higher centcr-of-niass energy, such as at the future e^e^ linear collider with ^/s = 500 
(LC500) GeV and 1000 (LCIOOO) GcV. We consider the case in which the lightest neutral Higgs 
boson of the NMSSM is composed dominantly of the singlet field. 

Let us first derive an upper bound on the mass of the lightest neutral Higgs boson in 
the NMSSM with explicit CP violation, where the Higgs singlet field contribute dominantly. 
The square of the radiatively corrected upper bound on the lightest neutral Higgs boson mass, 
"^fti, max i^'^ given as 

m| + (A^w^ - m|) sin^ 2/3 

A- ^2 

,2 ™2 



SttV (^2_^2^2 9{m^^, m-^J 




t2 tl' 

^ 3mi iXxcotpAi^+ At AjJ f m| ^ 3^4 
Att'^v'^ (m~^ — m~^) 

+ 87r2i;2 {ml -mlf ^^""^i' ""^^^ 

3mf (A,rtan.:?A; +.4, A;) / ™| \ 3/n? m| m? 
+ - . h log — + — ^ log(- 

(cos2/3A;^^ +sin2/3A;^2)2 [ m%mf.+ 



A.-^ A.1 \ A,l A.-^ / 



with 



Aj^ = At cos (f) + Xx cot P , 



12 



Aj^ = At + Xx cot (3 cos d> , 

Ag^ = cos (/> + Ax tan /? , (15) 
A^^ = At + Xx tan /? cos cj) . 

Note that A-^-^ and A^j already appear in the expression of 5Mf-. 

In the hmit of cos^ = cos(^c = 1) the above upper bound on the hghtest neutral Higgs 
boson mass reduces to the one obtained with CP conservation in the neutral Higgs sector. The 
first two terms appearing in the expression for m^^ are obtained from the tree level Higgs 
potential. They do not have any CP phase. The maximum value of A determines the tree level 
upper bound on the lightest neutral Higgs boson mass. The remaining terms in come 
from radiative corrections arising from the quark and scalar quark loops for the third generation, 
as well as the loops of W boson, the charged Higgs boson, and the charginos. 

In terms of m/^j max and m^j, it is useful to express the upper bound of the other neutral 
Higgs boson masses as [20] 



'^hz ^ '^hz, max 



m? 

hi, 


max 








1 - 






2 


max 


- (Ei=i Rl, 






1 - 








max 


- (Ei=i Rf, 






1 - 


(EliR^) 






max 








1 - 







™2 ^ 2 _ ,11, ma^ t/ ,n , 



"^/is ^ "^/is, max 



where Ri are defined as Ri = (On cos j3 + Oi2 sin/3) and satisfy the sum rule of Ei=i -^1 = 1- 

It is known that if the center of mass energy of the colliding electrons and positrons is larger 
than E'Y = mz + rrih^ (i = 1 to 5), the Higgsstrahlung process e^e~ Zhi (i = 1 to 5) is viable 
for the Higgs production. The cross section fij (i = 1 to 5) for hi production only through the 
Higgsstrahlung process is known to be related to that for the Higgs production in the standard 
model, (JSM, by 

(^j{mhj) = (^SM{mhj)Rj for j = 1 to 4 




(^5{mhs) = asuirrihs) \ 'i- - y Rj \ ■ (17) 



One can derive the parameter independent lower bound of cTj from the fact that cri(m/i^,max) ^ 
o-j (mhi)- 

In order to be systematic we assume that any of the five neutral Higgs bosons may be 
produced via the Higgsstrahlung process if kincmatically allowed. For each of them, we calculate 
the cross section for their productions. Let us dcnotcthe five cross sections as (Xj (i = 1 to 5); 
their dependence on the relevant parameters may collectively be expressed as 

(71 = ai{Ri,R2,R3,R4,mhi) , 

(Jfc = ak{Rl,R2,R3,R4,mhj,max) {k = 2to 5) . (18) 

By varying Rj (j = 1 to 4) and m/j^, one can obtain the maximum values for each cross sections, 
which would provide us a criterion on whether it is possible to detect any of them at the future 
e+e~ colliders. 
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For given Rj (j = 1 to 4) and m/i^, we may choose the largest one among the five cross 
sections: Let it be defined as 

(Jo = max{(7i, (72, (73, (74, (75} . (19) 

Now, we take the minimum of ctq by varying Rj {j = 1 to 4) and m/i^. This would tell us 
whether at least one of the five neutral Higgs bosons might be detected. Thus, we search ctq for 
^ Rj {j = 1 to 4) ^ 1 and for < m/j^ ^ ^hi, max) ^-nd set up its minimum value in order 
to examine the possibility of detecting at future e+e~ colliders at least one of the five neutral 
Higgs bosons of the NMSSM, at 1-loop level with explicit CP violation. 

In Fig. 3(a) we plot (Jq, as a function of Rj for ^/s = 500 (LC500) and 1000 (LCIOOO) GeV. 
For comparison, we also calculate the corresponding quantity for the case where the CP violation 
is absent. The solid curve is the result for the case of explicit CP violation, while the dashed 
curve is the result with no CP violation. Prom Fig. 3(a), one sees that for the case of explicit 
CP violation scenario, ctq can be as small as about 10 fb for ^/s = 500 GeV, and 2.5 fb for ^/s 
= 1000 GeV. On the other hand, for the case with no CP violation, the smallest value of ctq is 
about 16 for ^/s = 500 GeV and 3.0 fb for ^/s = 1000 GeV. By assuming an efficiency of 50%, 
the integrated luminosity at LC500 for 50 events would be required to be about 10 fb~^ in order 
to detect at least one of the five neutral Higgs bosons of the NMSSM at future e+e~ colliders if 
the CP symmetry is broken explicitly in the model. If there is no explicit CP violation in the 
NMSSM, the corresponding integrated luminosity needs to be 6.25 fb~^. 

In the literature, the worst case has been considered [21] when there is no explicit CP 
violation in the NMSSM: It is the case in which all the three neutral scalar Higgs bosons are 
degenerate in mass and Rl = Rl = Rl = 1/3. For our study, with explicit CP violation, the 
worst case would be i?? = 1/5 (i = 1 to 5). The five neutral Higgs bosons in our study may not 
be degenerate in mass. In Fig. 3(b), we plot a for i?? = 1/5 (i = 1 to 5). It is found that all of 
(7j (i = 2 to 5) except for ui increase as m^j increases. 

V. CONCLUSIONS 

We have investigated the neutral Higgs sector in the NMSSM with explicit CP violation at the 
1-loop level by using the effective potential method. For the radiatively corrected masses of the 
neutral Higgs bosons, we include not only the contributions of the quark and scalar quark for 
the third generation, but also the contributions of W boson, the charged Higgs boson, and the 
charginos. These additional contributions are found to be very crucial for the analysis of the 
neutral Higgs sector of the NMSSM, especially when the CP symmetry is broken explicitly. 

There are several complex phases in the NMSSM with explicit CP violation: At the NMSSM, 
three phases appear in the tree level Higgs potential and in the masses of the scalar top quarks 
and the scalar bottom quarks. We have assumed that these three phases are equal. At the 
1-loop level, besides these three phases, one additional CP phase arises from the 1-loop effective 
potential due to the contributions of the W boson, the charged Higgs boson, and the charginos. 

We have found that the mass of the lightest neutral Higgs boson may either increase by 16 
GeV or decrease by about 24 GeV, due to the contributions of W boson, the charged Higgs boson, 
and the charginos. This amount of variation in mass is quite significant to the relatively light 
rrih^ (25-40 GeV). Also the mixings among the neutral Higgs bosons are significantly affected 
both by the CP phase and the contributions arising from the contributions of the W boson, the 
charged Higgs boson, and the charginos. 
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The production cross sections for the neutral Higgs bosons in the NMSSM with expHcit CP 
violation are calculated via the Higgsstrahlung processes, e~^e~ — > Zhi (i = 1 to 5), when the 
lightest neutral Higgs boson is composed dominantly of the Higgs singlet field. The minimum of 
the production cross section for at least one of the five neutral Higgs bosons in the NMSSM is 
calculated to be about 10 and 2.5 fb for ^/s = 500 and 1000 GeV, respectively. These values are 
about one-third smaller than the corresponding numbers obtained by assuming CP conservation 
in the NMSSM Higgs sector. 
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FIGURE CAPTIONS 



FIG. 1. (a) : The plots of m^^, mh^, Im^J/mh^, and p as functions of the CP phase (pc- For 
m^^^ and rrih^, wc express them in GeV, and for |rn^_J/m/ij and p, in %. The vahies of other 
relevant parameters are = 7r/2, tan/3 = 3, A = A: = 0.5, = 800 GeV, Ak = At = 1000 GeV, 
X (= mr) = 500 GeV, mg = 900 GeV, and M2 = 300 GeV. 

FIG. 1. (b) : The same plots as in FIG. 1. (a), except for different parameter values = 7r/2, 
tan /? = 10, A = 0.03, = 60 GeV, A; = 0.4, Ak = 200 GeV, x = 500 GeV, mg (= mr = At) 
= 600 GeV, M2 = 100 GeV. 

FIG. 1. (c) : The same plots as in FIG. 1. (a) or (b), except for different parameter values = 
7r/2, tan/3 = 10, A = 0.7, ^a = 400 GeV, k = 0.5, Ak = 20 GeV, x = 50 GeV, mg = 500 GeV, 
mr (= At) = 1000 GeV, M2 = 100 GeV. 

FIG. 2. : The plot of rn^^ as a function of mj^^/irij^^, the ratio of the two chargino masses. The 
parameter values of the set of curves of feature (Al) are <p = it/2, tan l3 = 3, \ = k = O.b, A\ = 
800 GeV, Ak = At = 1000 GeV, x (= mr) = 500 GeV, and mg = 900 GeV. Those of feature 
(A2) are = tt/2, tan/3 = 10, A = 0.03, Ax = 60 GeV, k = 0.4, Ak = 200 GeV, x = 500 GeV, 
and mg = mx = At = 600 GeV, and those of feature (A3) are 4> = 7r/2, tan/3 = 10, A = 0.7, Ax 
= 400 GeV, k = 0.5, Ak = 20 GeV, x = 50 GeV, mg = 500 GeV, and mr = At = 1000 GeV. 
Each set consists of three curves corresponding to three different values of (j)c = (solid 
curve), 7r/2 (dashed curve), and tt (dotted curve). The extent of each curve corresponds to the 
variation of M2: We vary M2 from 1 to 500 GeV. 

FIG. 3. (a) : The minimum production cross section for at least one of the neutral Higgs boson, 
£70, is plotted as a function of Rj, for ^ = 500 (LC500) and 1000 (LCIOOO) GeV. The sohd 
curve corresponds to the case of the explicit CP violation while the dashed curve to no CP 
violation. 

FIG. 3. (b) : The production cross sections for each of the five neutral Higgs bosons in the 
NMSSM, (Tj (i=l to 5), are plotted as functions of m/^, for ^/s = 500 GeV. 
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FIG. 1. (a): The plots of m^_^, rrih-^, | J /m/i^ , and p as functions of the CP phase (pc- For 
m^^ and m/^j, we express them in GeV, and for |m^J/m/j^ and p, in %. The values of other 
relevant parameters are = 7r/2, tan/3 = 3, A = A; = 0.5, Ax = 800 GeV, Ak = At = 1000 GeV, 
X (= mr) = 500 GeV, mg = 900 GeV, and M2 = 300 GeV. 
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FIG. 1. (b): The same plots as in FIG. 1. (a), except for different parameter values = 
tan/3 = 10, A = 0.03, = 60 GeV, A; = 0.4, Ak = 200 GeV, x = 500 GeV, mg (= mr 
= 600 GeV, M2 = 100 GeV. 
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FIG. 1. (c): The same plots as in FIG. 1. (a) or (b), except for different parameter values = 
7r/2, tan/3 = 10, A = 0.7, = 400 GeV, k = 0.5, Ak = 20 GeV, x = 50 GeV, mg = 500 GeV, 
ruT (= At) = 1000 GeV, M2 = 100 GeV. 
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FIG. 2.: The plot of m^^ as a function of my^.^/mj^^, the ratio of the two chargino masses. 
The parameter values of the set of curves of feature (Al) are (j) = '7r/2, tan/? = 3, A = A; = 0.5, 
Ax = 800 GeV, Ak = At = 1000 GcV, x (= my) = 500 GeV, and mq = 900 GeV. Those of 
feature (A2) are (j) = 7r/2, tan/3 = 10, A = 0.03, = 60 GeV, A; = 0.4, Ak = 200 GeV, x = 
500 GeV, and mq = rriT = At = 600 GeV, and those of feature (A3) are <j) = 7r/2, tan/3 = 10, 
A = 0.7, Ax = 400 GcV, k = 0.5, Ak = 20 GeV, x = 50 GcV, mq = 500 GcV, and mr = At = 
1000 GeV. Each set consists of three curves corresponding to three different values of (j)c- 4>c = ^ 
(solid curve), tt/2 (dashed curve), and vr (dotted curve). The extent of each curve corresponds 
to the variation of M2: We vary M2 from 1 to 500 GeV. 
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FIG. 3. (a): The minimum production cross section for at least one of the neutral Higgs boson, 
ao, is plotted as a function of Rj, for = 500 (LC500) and 1000 (LCIOOO) GeV. The solid 
curve corresponds to the case of the explicit CP violation while the dashed curve to no CP 
violation. 
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FIG. 3. (b): The production cross sections for each of the five neutral Higgs bosons in the 
NMSSM, (Tj {i=l to 5), are plotted as functions of m/^, for -s/s = 500 GeV. 
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